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Abstract 

In this letter, we apply optimal control theory to design minimum-energy 7r/2 
and 7r pulses for the Bloch system in the presence of relaxation with constrained 
control amplitude. We consider a commonly encountered case in which the 
transverse relaxation rate is much larger than the longitudinal one so that the 
latter can be neglected. Using the Pontryagin's Maximum Principle, we de- 
rive optimal feedback laws which are characterized by the number of switches, 
depending on the control bound and the coordinates of the desired final state. 
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1. Introduction 



Optimal control theory [l| has been extensively used recently for the design 
of pulses that optimize the performance of various Nuclear Magnetic Resonance 
(NMR) and quantum systems limited by the presence of relaxation [1, IE IH H IE 
0, S H EE El El El, the dissipation due to random interactions between the 



system and its environment. In this letter, we employ tools from optimal control 
to derive minimum- energy 7r/2 and 7r pulses for a simple NMR system described 
by the Bloch equations. In particular, we study the case where transverse 
relaxation dominates the dissipation of the system and the control amplitude is 
bounded. 

The problem of optimal control of the Bloch equations and its closely re- 
lated corresponding problem for a two-level quantum system have received con- 
siderable attention. D' Alessandro and Dahleh [l4[ considered the problem of 
minimum-energy optimal control for a two-level quantum system without dis- 
sipation. Boscain and Mason [l5| examined the time minimal problem for a 
spin- 1/2 particle in a magnetic field neglecting dissipation. Sugny, Kontz and 
Jauslin Bonnard and Sugny Q and Bonnard, Chyba and Sugny [Ic| stud- 
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ied extensively the problem of time-optimal control for a dissipative two-level 
quantum system. 

In our recent work, we studied the problem of designing minimum-energy 
7r/2 and ir pulses for the Bloch system dominated by the transverse relaxation 



with unlimited control amplitude 12]. This dissipative system is of great prac- 



tical importance as it is a very good approximation for many applications of 
interest. In this article, we extend this previous work to consider the case where 
the control amplitude is limited, which accounts for realistic limitations of the 
experimental setup and also makes the problem more interesting from a control 
theoretic perspective. 

In the next section, we formulate the related optimal control problems of 
such pulse designs. The solutions of these problems are presented in Section 3, 
which is the main contribution of this article. Then in Section 4, we present 
some examples to demonstrate our analytical results. 



2. Optimal Control of Dissipative Bloch Systems 

In a resonant rotating frame, the Bloch equations with the longitudinal re- 
laxation neglected are of the form [l|| 



x = —Rx — UyZ 



-Ry + u x z, 



(1) 
(2) 
(3) 



where r = [x,y, z) is the magnetization vector, u x ,u y are the transverse com- 
ponents of the magnetic field and R > is the transverse relaxation rate. The 
above equations constitute a dissipative bilinear control system. By the follow- 



ing change of variables (see Fig. 1 (a) I and time rescaling 



a 

tan 9 
tan0 

tn,f=in 



= lnr = ln(y x 2 + y 2 + z 2 ) 



\fx?~-\ 

y/x 

Rtold, 



y 2 /z 



we arrive at a new system 



sin 2 9 



u — sm t) cos ( 
v cot 9 



(4) 
(5) 
(6) 



where u — (u x /R) sin</> — (u y /R) cos0, v = (u x /R) cos</> + (u y /R) sin</> are the 
normalized components of transverse magnetic field perpendicular and parallel 
to r± = (x,y), respectively. 

Note that v does not affect the angle 9 of the pulse. It just rotates r around 
z-axis, resulting in a waste of energy. Thus, optimality requires that v — and 
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(a) u _L rx (b) u || x-axis 



Figure 1: The optimal transverse magnetic field u is perpendicular to rx and its phase is 
constant (panel a). Without loss of generality, the experimental setup can be arranged such 
that u || x-axis. In this case, <j> = t/2 and r rotates in yz-plane (panel b). For convenience 
wc display z in the horizontal axis and y in the vertical, so that (r, 9) have the common 
configuration of polar coordinates on the plane. 



hence <fi = constant. If for example we choose u to be parallel to the a;-axis, 



then 4> = 7r/2 and r rotates in yz-plane, see Fig. 1(b) This is the case that 
we consider in this letter. Equations (0} and ([3]) are sufficient to describe this 
rotation. 

The optimal control problem that we would like to pursue is formulated as 
follows. Consider the dynamical system as in ([5]), starting from (r(0), 0(0)) = 
(1,0) (corresponding to (a(0),6>(0)) = (0,0)) and for a specified final value 
a(r) < a(0) = (equivalent to r(r) < r(0) = 1), what is the optimal bounded 
control u(t) with \u(t)\ < m E M + , < t < r, that accomplishes the transfer 
between the above starting point and the target point (a(r),0(r) = 7r/2or7r), 
while minimizing the energy E = u 2 {t)/2dtl Note that for the transfers that 
we study here, it must be m > 1/2, otherwise Eq. (|S|) reaches an equilibrium 
point #o < it/2. Also, the final time r is unspecified. 

3. Derivation of the Optimal Control 

The control Hamiltonian for the addressed problem is 

H = -u 2 /2 + \ e (u - sin 6» cos (9) - A Q sin 2 9 (7) 

where Xg, X a are the Lagrange multipliers. According to Pontryagin's maximum 
principle [l[, the necessary conditions for optimality of (it(£), 0(i), a(t), Xe{t), X a {t)) 
are 

X e = -8H/d9= A e cos 2(9 + A a sin 2(9 (8) 
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(a) Case of ir pulse (b) Case of tt/2 pulse 



Figure 2: The final point F[ can be reached by either following the minimum-energy path 
AF[, or traveling along AF\ up to point / and then leaving the system relax to F[ (panel 
a). Analogously, F^ can be reached by either following the minimum-energy path AF^, or 
traveling along AF2 and then leaving the system relax to FL (panel b). . 



A a = -dH/da = (9) 
u = argmaxi/ (it, 6, a, Xg, X a ) (10) 

u 

From (|9|), we immediately know that A a is a constant. Additionally, the optimal 
(u, 0, a, Xg, X a ) satisfies [l| 

H(u,e,a,Xg,X a ) = 0, 0<t<r. (11) 

Let E — mm„ JJ" u 2 (t)/2dt be the minimum cost corresponding to the op- 
timal solution. Using calculus of variations we find that small changes in the 
a-coordinate of the final point 5a T , and small changes in the final time St, 
produce the following change in the minimum cost 

SE = X a (r) 5a T - H{t) St. (12) 

Therefore 



A„(r) = dE/da T = dE/dr T ■ dr T /da T = dE/dr T ■ r>, (13) 

where r T — e clT is the radius of the final point. 
For the transfers that we examine here it is 



dE/dr T > 



(14) 



i.e. the larger the final r, the more energy is needed to quickly rotate the vector 
before it dissipates. To see this, we refer to Fig. [2j Let E, E' be the minimum 
energies necessary to reach the final points F\(r T , n), F[{r' T , tt), respectively, with 
The corresponding minimum-energy paths are AF% , AF[ . 



< ry, see Fig. 2(a) 



An alternative way to reach the point F{ is the following: travel along AF\ up 
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to the point /, with zj = Zp^ = —r' T , then set u = and wait until dissipation 
eliminates the y-coordinate, see (|3]). The energy E" spent for this travel is the 
portion of E necessary to reach /, so E" < E. By the definition of E' it is also 
E' < E" , and thus E' < E. Analogously, let E, E' be the minimum energies 
necessary to reach the final points F^ir-r , ir/2), -^(r^., 7r/2), respectively, with 



t' t < r T again, see Fig. 2(b) The corresponding minimum-energy paths are 
AF2,AF' 2 . An alternative way to reach the point F' 2 is the following: travel 
along AFi up to the point F%, then set u — and wait until dissipation brings 
the system at the point F%, see (J3|). The energy spent for this travel is the 
necessary energy to reach F2, i.e. E. Then, by definition of E', it is E' < E. 
Thus (TTJ} is true and from (TT5)) we have that A a (r) > 0. But A a = constant, so 
we can set 

A a = k 2 /2, k>0. (15) 

Having determined the sign of A Q , we first examine the case of unbounded 
control and then we use the developed intuition to study the general case of 
bounded control. 

3.1. Unbounded Control 

When the control u is unbounded, then from (|TU1) we conclude that dH/du = 
0. This condition and Eq. © yield 

u = \ e . (16) 
Using (fT6|) and (fT5|) . the condition (fTT|) becomes 

Ag — 2Ae sin 9 cos 9 — n 2 sin 2 9 — 0. (17) 
The optimal u is then given by the following feedback law 



u(9) = X e = sin6»(cos6»+ y/cos 2 9 + n 2 ). (18) 

Note that only the positive solution of the quadratic equation has physical 
meaning (corresponds to increasing 9) for the transfers that we study here. 
Using (fT6j) and (|T8l) . the validity of (jH) can be easily verified. Inserting (fT8| in 
([S]) we obtain the differential equation for the optimal trajectory 



9 = sin 9 a/cos 2 9 + k 2 (19) 
Eliminating time between Q and (fT9|) we obtain 

da sin 9 

d9~~ V^i 2 9 + K 2 ' 

Integrating the above equation from the starting point (0, 0) to the point (In r, 9) , 
we find the optimal trajectory 



... cos 9 + V cos 2 9 + k 2 . . 

ie) = i + vi + . 2 (21) 
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e(rad) ' ° "' " " " z 

(a) Optimal tt/2 pulse (b) Optimal trajectory 




e(rad) ' '* ' z 

(c) Optimal tt pulse (d) Optimal trajectory 



Figure 3: minimum-energy tt/2 (panel a) and tt (panel c) pulses for r T = 0.6. The correspond- 
ing trajectories are also shown (panels b,d). 



Setting 8 T = 7r/2,7r in the above equation, we find the optimal n for the tt/2 
and 7r pulses, as a function of the radius r T of the final point 



2r 

-tt/2 = -, ~2, «7r = : V ' ' (22) 

' 1 — rz 1 — r T 



The energy of the optimal pulses is calculated from 



Jo 2 J 29(6) 

using (JTHJ) and ([15]). The results for T = 7r/2,7r are 

E w/2 = -^, ^ = i±lL (24) 
1 — 1 — r T 

Using (f!M)) in (fT3")) . it is easy to verify the validity of (fT5|) with k given by (j2"2")l . 

In Fig. [3] we plot the optimal 7r/2 and tt pulses for r T = 0.6, as well as the 
corresponding trajectories. Observe that optimal u(9) is small close to 8 = 0, tt 
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(z-axis), directions that are protected against relaxation, while it is large close 
to 9 — tt/2 (y-axis), where dissipation is maximized and thus r must be rotated 
faster. 



3.2. Bounded Control 

We now move on to the case where the control is bounded, i.e. |w| < m, 
with to > 1/2 as pointed out before. The control Hamiltonian ([7]) is a quadratic 
form with respect to u that takes its maximum value at u = Xg, if \X$\ < to, and 
at the boundary point u — to if Xg > to. The other boundary point, u = —m, 
corresponds to decreasing 8 and has no physical meaning for the transfers that 
we examine here. Initially, the situation is as in the previous case where the 
optimality condition u = Xg holds and the optimal control is given by (I18[) . 
Angle 9 increases following (TTO1) and u, Xg change accordingly. Now suppose 
that at some point the control reaches the maximum allowable value to. From 
(fl"8|) we see that this happens at the angles that satisfy the equation 



sin#(cos0 + v cos 2 9 + k 2 ) = to, (25) 

which is equivalent to the following quadratic equation for cot 9 

2 k 2 

cot 2 9 cot 6 + 1 ~ = 0. (26) 

to m 2 

If #i,#2 are the solutions of (|2"6"|) in [0, 7r], then it is easy to show that for 
9 £ (9i,92) the following inequality holds 

sin 0(cos 9 + \/cos 2 9 + k 2 ) > to. (27) 

In the interval 9 G (^1,^2) where the above inequality is true, the relation 
u = Xg gives u{9) > m which is not permissible. Thus, the optimal control in 
this interval is u(9) — to. From (|lip we can find the Lagrange multiplier Xg for 
the same interval, which is 

_ to 2 + k 2 sin 2 9 
Xe{9) ~ 2(m- sin0 cos 9) (28) 

and Xg{9\) = to from (|16p . This Xg satisfies the optimality condition ([8j with 
9 evolving in time according to ([5]) and u(0) = to. It is not hard to verify using 
l[2"7J). ([25|1 that Xg > to for 9 e (9i,0 2 ), thus u{9) = m is indeed the optimal 
control in this interval. We call the solutions 81,62 of (|26|). where a change in 
the optimal control occurs, the switching angles. After the second switching, 
the optimal control is given again by (|18|) . with the same K as in the initial 
phase, since A a = k 2 /2 is constant along the optimal trajectory. 

As explained above, the switching angles determine the optimal feedback 
law. For to > 1, Eq. (|2"6"|) has real solutions for k > \Jm 2 — 1. For 1/2 < to < 1 
it has real solutions for every n > 0. We examine separately these two cases 
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Figure 4: Switching curves AB and BC\ for the case m > 1. The outermost curve AC\ 
defines the reachable set from the starting point ^4(1,0) when u < m. Curve BC2 is the 
optimal trajectory for k = \/m 2 — 1 and 9 £ [6b, 7r]. The perpendicular axis 6 = tt/2 crosses 
the curves ACi, BC\, BC2 at the points D\, D2, D3, respectively. 



3.2.1. m > 1 



For k > \/m 2 — 1 there are two switching angles, given by 



, / 1 ± Vre 2 - m 2 + 1 \ 
^ 2 = cot-^ j (29) 

In this letter the range of the function cot -1 is considered to be [0,7r], so 

cot _1 (a;) = 7T - cot-^-x), x < 0. (30) 



In the case k = \/m 2 — 1, the two angles obtain the common value Ob = 
cot _1 (l/m). The angle of the first switching point is in the range 6\ £ [0,9b], 
while that of the second 62 E [Ob , 7r] . We find the equation of the first switching 
curve, i.e. the curve composed by the points (rx,#i). Before the switching the 
optimal control is given by (|18l) . so each of the points (ri,#i) belongs to an 
optimal curve of the form (|2"Tj) . i.e. 



cos 6»i + Vcos 2 6»i + k 2 

n(di) = r— -f • (31) 

1 + vl + re 



But this k is related to the switching angle 0\ through (1291) . which we re- write 

as 

k 2 = (m cot 0i - l) 2 + m 2 - 1 (32) 

since only k 2 appears in pip . These two equations determine the first switching 
curve, with the angle 61 in the interval 9\ € [0, 63]- 

After the first switching, the optimal control takes the value u(9) = m. It 
maintains this value until the second switching, at angle 62- Between the two 
switchings the evolution is given by 

da sin 2 9 , , 
= (33) 

d9 m — sin 9 cos 9 ' 
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as we derive from ([5]) with u — m. Integrating the above equation from 9\ 
to 62 we find the radius of the second switching point 



ra(0 2 ) = ri( 



/2m- sin 20i 
7l) ^ 2m-sin2fl 2 6XP 



/(fi,0 2 ) 



where 



and 



f(9 1 ,0 2 ) = cot- 



( 2m cot 6» 2 -l 
V V2" 



1 



— cot 



V4m 2 - 1 



2m cot #1 — 1 



V47 



1 



cot _1 (2/m - cot 6»i). 



(34) 



(35) 



(36) 



Therefore, to every first switching point (ri,0i) corresponds a second switching 
point (r 27 9 2 ) with angle given by ([31)1) an d radius given by ([34"]) . ([33)1 . These 
points compose the second switching curve. The two switching curves are plot- 
ted in Fig. HI curves AB and BC\. The joint point is B(tb,9b), where 



r B = 



V 'm 2 



1 



m + 1 

The outermost curve ^4Ci corresponds to the trajectory traveled for u(9) 
m,9 £ [0, 71"]. The equation for this trajectory can be found by setting (ri, 9\) 
(1,0) (starting point A) in (J33]), ([55]). It is 



(37) 



r 3 (0 3 ) 



2m 



2m — sin 2#- 



■ exp 



1 



V47 



: COt 



2m cot 9"i 



V47 



(38) 



where we used (ra, 9^,) to denote a point on the curve and 6*3 € [0, 7r]. The points 
between this curve and the horizontal axis define the reachable set from A for 
a specific control bound m. This curve meets the axes 9 — tt,tt/2 at the points 
Ci(rc 1 ,7r),Di(r Dl ,7r/2), where 



r D 1 



exp 
exp 



V4T 



V4j 



1 



7T — COt 



V47 



1 



(39) 
(40) 



These are the points with the largest radius along these axes, that can be reached 
from A(1,0). Using ([31]), ([32]), (j3J}, ([35]) and ([M]) we find that the second 
switching curve BC\ crosses the axis 9 = tt/2 at the point D2(td 2 , where 



rD 2 



y/m 2 + 2 



1 + Vm 2 + 1 



exp 



1 



V4w? 



: COt 



1 



V4m 2 



(41) 



As we mentioned above, switching takes place only for n > Vm 2 — 1. For 



1 there is no switching and the optimal trajectory is given by (|21[) . 



9 



For these values of k, the optimal trajectory crosses the axis 9 — n, tt/2 at points 
with radius 

Vk 2 + 1 - 1 K . . 

r7r = / 9 ; 7' r V 2 = / 9 = 7 ( 42 

Vk 2 + 1 + 1 \/k 2 + 1 + 1 

respectively, both increasing functions of k > 0. In Fig. 2] we plot the optimal 
trajectory without switching for the largest permissible value k = V m 2 — 1 



cos 6» + Vcos 2 O + m 2 - 1 

r o( d o) = —j (43 

m + 1 

and for 9o £ [#b, 7r]. It crosses the axes 9 = tt, tt/2 at the points Ci(jc^ , tt), D^{td 3 , 
where 

m — 1 v m 2 — 1 , 

r C 2 = — —r, r Da = — — 44 

m + 1 ?7i+l 

These are the largest radius points along these axes, that can be reached without 
switching. 

Using the construction shown in Fig. [4] we can find the switching points and 
the optimal control for any final point of the form Fx(r T ,ir) or ^(rv, 7r/2). If 
Fi £ C1C2 then the optimal trajectory after the second switching is 



ia\ ra \ cos 9 + Vcos 2 9 + K 2 
r(9) = r 2 {9 2 ) — — == (45) 

cos 9 2 + V COS^ 82 + K 

where 

k 2 = (1 _ mcot6» 2 ) 2 + m 2 - 1 (46) 

and S 2 (r2, 62) is the second switching point. Eq. (|4"5j) is found after integrating 
(HOD from 5 2 to Fj, while (@BJ holds because S2 € BCi, see flU) . The final 
point Fi(r T ,x) belongs to this curve, so we find 



cos 6» 2 + V cos 2 2 + k 2 

r2(gzHr - -i + vi + . 2 (47) 

Plotting this curve for 9 2 £ [6b, 7r], it crosses the second switching curve at the 
second switching point S2. The first switching angle can be found from 

9 1 =cot- 1 {2/m-cot9 2 ). (48) 

Having determined the two switching angles and the optimal k, the optimal 
feedback control is also determined. If F\ £ C2O then no switching is necessary 
and the optimal control is given by (ITS)) with n — K n given by (l22t . 

For F 2 £ DiD 2 there is only one switching. The optimal trajectory after the 
switching is 



2m- sin26»i 
r(g)=n(*iW 2m _ sin2g ex P 



\/Am 2 - 1 



(49) 
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Figure 5: Switching curves AB\ and B2C1 for the case 1/2 < m < 1. The outermost curve 
ACi defines the reachable set from the starting point A(l, 0) when u < m. Curve B2O is the 
optimal trajectory, after the second switching, for the lowest value ft = 0. The perpendicular 
axis 6 = 7r/2 crosses the curves ACi, B2C1 at the points D\,D2, respectively. Note that for 
these values of m, switching is inevitable for the transfers that we examine. 



The final point F2(r T ,ir /2) belongs to this curve, so we find 

'M, tt/2) 



n(0i) 



2 m 

2m — sin 2^i 



exp 



%/4to 2 - 1 



(50) 



Plotting this curve for 6\ £ [0,6b], it crosses the first switching curve at the 
switching point S\(ri, 9\), so we find the switching angle and the corresponding 
optimal k. If Fi £ -D2-D3 then there are two switchings and the situation is 
similar to the one previously described, with 



»~2(02 



COS 6 2 + %/ COS 2 6 2 + K 2 



(51) 



instead of (|47p. since = 7r/2 at the final point. If F2 £ D3O then no switching 
is necessary; the optimal control is given by (|18l) with k = k^/2 given by (|22|) . 



1/2 < m < 1 

For this case the two switching curves are described by the same equations 
as before but now there is no common point but a gap between them, due to 
the fact that the two roots of (f26| are distinct for all n > 0, see Fig [5] The 
range for the switching angles is 9\ £ [0,0 3^,62 £ [0b 2 j 7F ] where 



cot 



cot 



1 - VT^ 



(52) 



(53) 



The switching points B\ , B2 correspond to the case k = and their radius is 

r Bi = \ o ( 54 ) 
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tb 2 = r Bl x exp 



: COt" 1 



V4m 2 - 1 \\/4m 2 - l-y/T 



2tti 2 



(55) 



For the lowest value re = we plot the optimal curve after the second switching. 
It is 

fa \ cos9 ° fKR\ 
ro{9 ) = r B2 j. — (56) 

cos 0b 2 

This curve meets the axis 9 = ir/2 at the origin O. Points C\,D\, Z?2 are defined 
as before and their coordinates are the same. The situation is as depicted in 
Fig. [51 Observe that for 1/2 < m < 1, switching is inevitable for the transfers 
that we examine; since the control is more restricted than the previous case, 
the boundary value has to be used to achieve the desired transfers. Any final 
point of the form Fi(r T ,ir) belongs to C\0 and there are two switchings in the 
corresponding optimal trajectory. The same happens for i<2(r T ,7r/2) G D2O, 
while for F2 G -D1-D2 there is only one switching. 

3.2.3. Summary of the results 

We summarize the above results 

• to > 1 

1. 9 = 7T 

(a) rc 2 < r T < rc x , two switchings 

(b) r T < rc 2 , no switching 

2. 9 = ir/2 

(a) ru 2 < r T < rr> l , one switching 

(b) 7\d 3 < r T < rr> 2 , two switchings 

(c) r T < ro 3 , no switching 

• 1/2 < 771 < 1 

1. 9 = 7T 

(a) r T < ret , two switchings 

2. 9 = tt/2 

(a) ru 2 < r T < rr> l , one switching 

(b) r T < ro 2 , two switchings 

If 0i, 6>2 are the switching angles then the optimal control for 9 G (0i, #2) is 

u(9) = 771 (57) 

while outside this interval is 



u{9) = sin 9 (cos 9 + y/cos 2 9 + k 2 ) (58) 

where 



K = \J{m COt 6>1 - l) 2 + 771 2 - 1 = - 77lCOt6> 2 ) 2 + m 2 - 1. (59) 

The switching angles are calculated as described above and are related through 

cot 6»i +cot6> 2 = 2/m. (60) 
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4. Examples 



Here we present some examples using specific values for the upper bound 
m and the coordinates (r T ,0 T ) of the final point. We start with the case m = 
2,r T = 0.39, 9 T = it. Using the results of the previous section, it is not hard 
to see that the optimal trajectory contains two switching points. Plotting Eq. 
(|4"T|) we find that it crosses the second switching curve at the switching point 
S2 with angle 9 2 = 1.7766 rad. The first switching angle is found from (IM|) to 
be 6»i = 0.6912 rad. The optimal cont rol is given by ([57) , {5SJ with « = 2.2382, 
as it is determined from (1B"9|) . In Fig. 6(a)| |6(b"J| we plot the optimal pulse and 
the corresponding optimal trajectory. The two switching points Si , S2 are also 
shown. 

For the next example we use the values m = 2, r T = 0.61, 9 T — ir/2. In this 
case, the optimal trajectory contains only one switching point. Plotting Eq. (|50[) 
we find that it crosses the first switching curve at the switching point Si with 
angle 0i = 0.6124 rad. The optimal control is given by (l58t with k = 2.5322 for 
9 e [OA] and by J57J) for 9 e (0i,tt/2]. In Fig. [6(c)| [6(d)1 we plot the optimal 
control and the corresponding optimal trajectory, with the one switching point 
Si. 

The last case that we consider is m = 0.95, r T = 0.2, 9 T — ir/2. The optimal 
trajectory contains two switching points. Plotting Eq. (|5ip we find that it 
crosses the second switching curve at the switching point £2 with angle 62 = 
1.1456 rad. The first switching angle is found from (|60l) to be 6\ — 0.5442 rad. 
The optimal control is given by ([ST|) . (jS"5|) with k = 0.4766, as it is determined 
from ([Sn]). In Fig. 6(e)| |6(f)| we plot the optimal pulse and the corresponding 
optimal trajectory. The two switching points Si,S2 are also shown. 



5. Conclusion 

To conclude, in this letter we calculated minimum- energy tt/2 and ir pulses 
for Bloch equations in the case where transverse relaxation dominates and the 
control amplitude is bounded, using optimal control theory. This work is ex- 
pected to find applications in NMR Spectroscopy, Magnetic Resonance Imaging 
(MRI) and Quantum Information Processing, serving as a reference for numer- 
ical studies of more complicated and realistic situations that incorporate for 
example longitudinal relaxation and magnetic field inhomogencity. 
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Figure 6: minimum-energy pulses (a, c, e) and corresponding optimal trajectories (b, d, f) for 
various values of the control upper bound m and the coordinates (r T ,8 T ) of the final point. 
The switching curves and points are also shown. 
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